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Abstract 

Evidence of a primitive spinodal decomposition has been obtained for 
central Ni+Ni Heavy Ion Collision, since higher order charge correlations 
show a peak when four fragments of size equal to 6 are produced with 
an excitation of 4.75 MeV. This can be considered as a signature of a 
primitive breakup in equal sized fragments with a privileged fragment size. 
This computational result confirms other experimental and theoretical 
evidences about spinodal decompostion in HIC. 


1 Introduction 

Heavy Ion Collisions are expected to experience a liquid-gas phase transi¬ 
tion due to the specifics of the internucleonic interaction, which is attrac¬ 
tive in both the long and intermediate ranges and repulsive in the short 
range. It is possible that a wide zone of phase space is explored when 
two nuclei collide, including a region where liquid and gas phases coex¬ 
ist, namely the spinodal region. At this spinodal region, incompressibility 
is negative and uniform nuclear matter is unstable, leading to multifrag¬ 
mentation due to the increase of density fluctuations D 0 Dynamical 
simulations based on Boltzmann equation, such as Landau-Vlasov (LV), 
Boltzmann-Uehling-Uhlenbeck (BUU) or Boltzmann-Nordheim-Vlasov (BNV), 
describe the time evolution of the density of a one-body system, ignoring 
those correlations whose order is larger than the order of binary correla¬ 
tions, neglecting fluctuations around the mean trajectory of the system, 
which altogether comes out to be quite inconvenient to study the spinodal 
instability zone [T. As shown by Guarnera et. al. [[]], a spinodal de¬ 
composition produces a ’’primitive breakup” where equal sized fragments 
have a privileged fragment size, which is related to the wave lenghts of 
the most unstables modes of nuclear matter. Tabacaru et. al. have ob¬ 
tained reduced velocity correlations between fragments and did not find 
a bubble-like profile, which excludes surface instabilities that might cause 
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multifragmentation. They also computed higher order charge correlations, 
defined as the ratio between the number of correlated fragments and the 
number of uncorrelated fragments, obtaining evidence of a privileged pro¬ 
duction of equal sized fragments |1J. Chomaz et. al. have introduced a 
scenario inspired on experimental data, where a gently compressed sys¬ 
tems expands and reaches thermal equilibrium approximately at the time 
when the system enters into the spinodal region. At this moment, density 
fluctuations break up the system into several hot fragments and parti¬ 
cles. When fragments are released from the nuclear force, configuration 
frozens and fragments only interact with each other via coulomb force. 
At this moment, system has explored so much phase space that it can be 
described by statistical models and there is no contradiction between sta¬ 
tistical or dynamical approximations. Statistical approximations describe 
the time evolution of the system and the phase diagram. Meanwhile, 
dynamical approximations start up in phase diagram and are more re¬ 
lated to the thermodynamics of non extensive systems. Barranon et. al. 
have obtained computational evidence about the inverse relation between 
entropy and the residual size, using LATINO dynamical model to study 
the spinodal decomposition region of central Ni+Ni HIC at intermediate 
energies [5] . In this very study, evidence is obtained about a primitive 
spinodal decomposition for Ni+Ni central HIC at intermediate energies. 


2 Methodology. 

Heavy Ion Collisions were simulated using LATINO semiclassical model 
where binary interaction 0 is reproduced with a Pandharipande potential 
built up of Yukawa potentials linear combinations, whose coefficients are 
designed to both reproduce nuclear matter ground state properties and 
to fulfill Pauli exclusion principle J7],: 


Vnn = V pp = Vo(- 


(!) 


and: 

/e -Mr7 ' e _Mrrc x , e _ ''“ r e _M “ r ' a x 
V np = V r (— -— )-Va(— -—) (2) 

Clusters are identified with an Early Cluster Recognition Algorithm 
that optimizes configurations in energy space. A most bound partition is 
obtained minimizing the sum of energies of the clusters belonging to each 
partition: 


{Ci} = argmin [E {Ci} = E E% t ] (3) 

i 

where the energy of each cluser is given by : 

E %t=Y,iY, K ? M+ E M ( 4 ) 

i ij€.Ci j,k£Ci,j<k 

where the first sum includes the partition clusters, Kf M is the kinetic 
energy of the particle j measured in the center of mass of the cluster 
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containing particle j, and Vij is the internucleonic potential. The algo¬ 
rithm uses ’’simmulated annealing” to find the most bound partition and 
optimizes the partition in energy space. 

Projectile energy is in the range of 600 to 2000 MeV and system evolves 
until its microscopic composition is rather frozen though some monomers 
are ejected. This time can be identified using the Persistence Microscopic 
Coefficient, defined as the probability that two particles belonging to the 
partition X remain bound in partition Y: 


P[X,Y) 


1 

cluster 



cluster 


(5) 


where bi is equal to the number of pairs of particles belonging to the 
cluster Ci of partition X while a; is equal to the number of particle pairs 
belonging to cluster Ci of partition X that also belong to a given cluster 
C'i of partition Y. m is the number of particles in cluster Ci. 

Higher order charge correlations were introduced by Moretto et. al. 
f5; and are given by: 


Y(AZ,< Z >) 

Y’ (AZ, <Z>) 

where Y^AZ, < Z >) is equal to the number of fragments produced for 
given values of A Z and < Z >. 


( 6 ) 
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3 Results 

Higher order charge correlations were obtained for central HIC Ni+Ni. 
These higher order correlations show a peak for a fragment size equal 
to 6, with four equal sized fragments produced with this privileged size 
and an excitation equal to 4.75 MeV. This can be considered a primitive 
breakup signature in equal sized fragments of a privileged size. Hence, this 
dynamical study supports other experimental and theoretical evidences of 
spinodal decomposition in HIC, such as the experimental study reported 
by Borderie et. al. [SJ where a fosil spinodal decomposition signature for 
HIC collision 129 Xe -\- rlat Sn at 32 AMeV was obtained. At that very ex¬ 
periment, liquid-gas coexistence was observed and evidence was obtained 
of a first order phase transition for a non extensive particle system. This 
HIC phase transition has been identified in previous studies with several 
signatures, namely Campi scattered plotts m as well as fragment size 
distribution power laws El. 

4 Conclusions. 

Once LATINO dynamical model was applied, higher order charge correla¬ 
tions provided computational evidence about the spinodal decomposition 
in the early stage of fragmentation for central HIC Ni+Ni at interme¬ 
diate energies. This confirms previous evidence reported by others of a 
fosil spinodal decomposition, indicating liquid-gas coexistence and a first 
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